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ABSTRACT 


CURRENTS BETWEEN TETHERED ELECTRODES IN A 
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Laboratory experiments on important plasma physics issues of 
electrodynamic tethers are performed. These include current propagation, 
formation of wave wings, limits of current collection, nonlinear effects and 
instabilities, charging phenomena, and characteristics of transmission lines 
in plasmas. The experiments are conducted in a large afterglow plasma (100 
cm x 200 cm, n < 10 12 cm* 8 , kT < 3 eV, B < 100 G, Ar, 3 x 10*“ Torr) . The 
current system is established with a small electron-emitting hot cathode 
tethered to an electron-collecting anode, both movable across the magnetic 
field and energized by potential differences up to V « 100 kT e /e. The total 
current density in space and time is obtained from complete measurements of 
the perturbed magnetic field, J — V^x B(r,t)/#i . The fast spacecraft motion 
is reproduced in the laboratory by moving the tethered electrodes in small 
increments, applying delayed current pulses, and reconstructing the net field 
by a linear superposition of locally emitted wavelets. With this technique, 
the small -amplitude dc current pattern is shown to form "whistler wings" at 
each electrode instead of the generally accepted "Alfven wings." For the beam 
electrode, the whistler wing separates from the field-aligned beam which 
carries no net current. Large amplitude return currents to a stationary 
anode generate current -driven microinstabilities, parallel electric fields, 
ion depletions, current disruptions and time-varying electrode charging. At 
appropriately high potentials and neutral densities, excess neutrals are 
ionized near the anode. The anode sheath emits high-frequency electron 


(NASA-CR-186128) CURRENTS BETWEEN TtTHERED 

electrodes in a magnetized laboratory plasma 

(California Univ.) 70 p CSCL 201 


N91- 12 344 


Unci as 

G3/75 0252156 



transit- tine oscillations at the sheath-plasma resonance. The beam generates 
Langmuir turbulence, ion sound turbulence, electron heating, space charge 
fields, and Hall currents. An insulated, perfectly conducting transmission 
line embedded in the plasma becomes lossy due to excitation of whistler waves 
and magnetic field diffusion effects. The implications of the laboratory 
observations on electrodynamic tethers in space are discussed. 

1 . Introduction 

Electrodynamics is of general interest in plasma physics [Shercliff, 

1965] and of particular importance in space plasmas when conducting objects 
are in relative motion with respect to a magnetized plasma [Drell et al., 
1965; Goldreich and Lynden-Bell, 1969]. The possibility of tethering 
satellites with long wires has stimulated many proposed active experiments 
[tfilliamson et al., 1978; Gross! and Colombo, 1978; Dobrowolny, 1978). If a 
tether wire of length and direction T moves with velocity v across the 
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earth's magnetic field B q a motional EMF, V - (v x B q ) ■ i, is Induced which 
could drive a current In a circuit involving the wire, the end electrodes and 
the ambient stationary plasma. Such a system could be used to generate power 
and thrust [Banks et al., 1980; Martinez -Sanchez and Hastings, 1987], inject 
electron beams [Winckler, 1980; Obayashi et al., 1984], generate low 
frequency waves [Belcastro et al., 1982; Rasmussen et al., 1985; Barnett and 
Olbert, 1986; Hastings et al., 1988], and perform remote sensing [Colombo et 
al., 1975]. Important issues for the envisioned uses of tethers in space are 
its mechanical properties such as drag, insulation, wire breakage, etc., and 
the electrical behavior which includes limits on current collection at the 
electrode -plasma interface (contactors), current propagation path ("Alfven 
wings") and current closure across B^, Ohmic and radiation resistance of 
tether and contactors, and nonlinear plasma phenomena (instabilities, double 


- 2 - 


layers , ionization) . 


While most of these problems have been addressed in theoretical studies, 
there are likely to be many unforeseen or unknown effects which can only be 
found from actual experiments. Some experiments on short tethers have been 
performed with rockets [Williamson et al., 1981, Myers, et al., 1989] and 
selected topics such as contactors are under investigation in the laboratory 
[Wilbur and Laupa, 1988], but the major electrodynamic issues of tethers 
remain untested. In particular, considerable attention has been given to 
the type of electromagnetic wave that a tether and its contacting ends 
might emit. While most of the works have concentrated on the Alfven regime, 
as first done by Drell et al., [1965], there has been some inference that 
emission will also take place at other frequencies. Rasmussen et al., [1985], 
for example, speculated that shear-Alfven waves and whistlers will also be 
emitted leading to energy loss on a cone centered on the magnetic field lines 
intersected by the tether's contacts. Barnett and Olbert [1986] and, more 
recently, Hastings et al., [1988] have shown that wave emission in the bands 

< « < O e and «p < w < wuh can take place and, with appropiate 
parameters, dominate that of the Alfven band (w < Oj). 

Recent laboratory experiments on magnetic reconnection [Stenzel et al., 
1986] have shown the feasibility of detailed observations of dynamic current 
systems in plasmas. In the present experiment we investigate a current 
system set up with tethered electrodes. We show how the current is 
distributed in space and time [Urrutia and Stenzel, 1986, 1989] resulting in 
the formation of "whistler wings" [Stenzel and Urrutia, 1989], describe 
instabilities which arise due to inherently nonlinear circuit elements 
[Stenzel and Urrutia, 1986; Stenzel, 1988], and discuss questions of scaling, 
boundaries and collisions. While the laboratory experiment cannot model all 
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aspects of tethers In space (e.g. the MHD far-zone or the motional EMF) It 
does support certain predictions (e.g. wave wings), contradicts others (e.g. 
distant current closure), and reveals many interesting new phenomena (e.g. 
current disruptions, current-neutralized beams, wing spread). 

After the description of the experimental setup and measurement 
techniques we first present the results of small amplitude currents which 
allow us to construct whistler wave wings based on linear superposition 
principles. Next we show the nonlinear effects and instabilities, and last 
discuss the implications of the laboratory results on tethers in space. 


2 . Experimental Setup and Measurement Techniques 

The experiment shown schematically in Figure 1 is performed in a large 
volume (1 m diam. , 2 m length) uniform (n/Vn “ 10 n), quiescent (£n/n m 1%), 
weakly collisional (y^/ci^ *10 ) afterglow plasma of density n^ <10 cm , 

temperature kT^ ** 10 kT^ < 3 eV, produced by a repetitively pulsed dc 
discharge (40 V, 800 A, t^ « 3 m, t * 0.5 s) in Argon at a pressure 
p “ 3 x 10 4 Torr. The quality of the plasma is due to several factors, (i) 
the uniform plasma production with aim diam. oxide-coated cathode [Stenzel 
and Daley, 1980], (ii) the uniform dc axial magnetic field (B q < 100 G) , and 
(iii) the simple linear geometry. The Maxwellian afterglow plasma is best 
suited for studying imposed current systems since the medium is inherently 
current-free , and for understanding nonlinear dynamic processes since the 
plasma has no volume sources or sinks. The density decay due to surface 
recombination in a large chamber is slow (t^ r 1000 ps) compared to the 
time scale of the current propagation and closure (t <: 5 /is). 

In analogy with a tether system in space, we insert into the plasma an 
electron emitter (oxide cathode, 1...3 cm diam.) and an electron collector 
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(metal disc or sphere, 1...3 cm diam.), well separated across (dr 15 cm) 

electron Larmor radius, r. <1 cm) and both electrodes movable together 

Le 

across i . However, in contrast to the tether system in space, the potential 
o 

difference between the electrodes is not generated by a motional EMF but by 
an external pulsed power supply. In general, the electrodes are electrically 
floating with respect to the chamber wall, hence find their self-consistent 
potential relative to the plasma potential as in space. However, in the 
laboratory we have the option to close the current circuit from either 
electrodes to the end chamber wall (equivalent to the collisional ionosphere) 
which is very helpful in isolating nonlinear effects of the floating current 
system. 

Different diagnostic tools are used to measure fields and particles. A 
vector magnetic probe consisting of three identical, orthogonal, 
electrostatically shielded loops (1 cm diara. , 30 turns) is used to measure 
the magnetic field perturbations caused by the plasma current. The probe can 
be moved In three dimensions (x,y,z) and resolves small fields (A B * 10 8 G) 
with high time resolution (At ** 0.1 p s). A complete field measurement, 
B*(x,y,z,t), typically involves three simultaneous recordings of the vector 
components versus time, ensemble averaged over repeated events (N > 10) so as 
to Improve the signal -to -noise ratio. The spatial position is then varied in 
the plasma volume containing the current system (up to 5000 locations). The 
digitized data is stored and subsequently evaluated to yield the complete 
vector field of the total current density j*(x, y,z,t) - V x B/a* 0 , without any 
assumptions about field symmetries, and verifying V • B * 0 to check the data 
quality. Surface integrals of the current density yield net currrents which 
can be compared with the measured tether current so as to estimate where the 
effective current closure takes place. We note that the total current 
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density 7x1? includes the displacement current 8D/3t, and is always 
divergence * free (7 • J “0). Therefore, field data alone do not identify the 
transport mechanism. 

Plasma density, temperature and potential are obtained from small movable 
Langmuir probes. Due to the rapid changes in plasma parameters, the full 
probe characteristics are not obtained in the traditional way (full sweep 
within a short time). Instead, the bias of the probe is stepped through the 
entire range over succesive shots as documented in Urrutia and Stenzel [1986] 
and, more recently, by Sheridan and Hayes [1988]. Electrostatic fluctuations 
are detected with coax- fed rf probes and analyzed by cross -correlation 
techniques [Urrutia, 1987]. High-frequency instabilities are observed with 
microwave antennas at the plasma boundary and within the plasma volume 
[Stenzel, 1989]. The electron beam (100... 500 eV, 1 — 10 A) generates 
visible light by exciting atomic transitions in ions and neutrals [Griem, 
1964]. The beam trajectory is externally observable since the cold 
background plasma is dark and optically thin. The space -time dependence of 
pulsed electron beams has been measured using a fast (10 ns) photomultiplier 
detector with collimated optics, [Chan and Stenzel, 1988]. The propagation 
speed of particles and currents is, in general, different. 

3. Linear Currents 

It is appropriate to start with small -amplitude currents which do not 

give rise to instabilities or plasma modifications. Typically, this holds for 

current densities on the order of the ion saturation current densities, J - 

ne(kT /m,) 1 / 2 We are interested in the propagation speed for time-varying 
e i ■ 

currents, the propagation path of currents between stationary and moving 
electrodes, and the plasma currents induced by an insulated tether wire. 

3.1 Current Propagation 
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Although in the frame of a uniformly moving tether the induced voltage 

and current are time independent, the stationary plasma experiences a current 

pulse on each flux tube due to the finite time an electrode stays on a field 

line. Thus, the plasma currents are inherently ac currents, and the basic 

problem is how a current pulse propagates away from an electrode. 

Figure 2 shows the time dependence of the perturbed magnetic field 

B,(±B* ) at different axial distances z (IB* ) from a disc electrode which 
^ o o 

collects electrons for a pulse of At * 30 ps duration. The magnetic field 
exhibits an oscillatory transient, both at tum-on and turn-off, with 
frequency decreasing in time and propagation speed decreasing with frequency. 
The typical propagation speed (1.25... 5 x 10 ? cm/s) of the magnetic field B^, 
hence axial current density J - (^ r) 1 3(rB,)/3r, is smaller than the 

Z o P 

electron thermal speed (v - 6 x 10 T cm/s, kT - 1 eV) . Furthermore, the 

6 C 

propagation speed is observed to decrease with higher electron densities n^ 
and lower dc magnetic fields B^, while being essentially independent of the 
electron temperature. Thus, the current does not propagate at the random 
speed of the collected particles but is carried by a wave which is identified 
in the following experiment. 

In order to excite a single mode rather than a broad wave packet, the 
current pulse is shortened and repeated continuously (e.g., f - 1 MHz, square 
wave). Via an interferometer circuit, the fundamental 1 MHz magnetic field 
perturbation is selected and displayed vs. axial position in Figure 3a. A 
variable time delay r is inserted into the signal path in order to determine 
both magnitude and sign of the phase velocity v ^ " w/k - -z/r from the 
interferometer output, B^cosCwr + kz). The results of Figure 3a show an 
undamped magnetic wave propagating along B^ away from the electrode located 
at z - 0. The wave polarization is obtained by recording two orthogonal 
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field components (B ,B ) vs. z, as shown in Figure 3b, and constructing the 

x y 

vector £ at different positions z (see insert). The wave is found to be 

right-hand circularly polarized. Finally, the dispersion curve u> vs. k| is 

obtained by measuring wavelengths at different frequencies, displayed in 

Figure 3c in normalized form. Within measurement accuracy, the data points 

fall on the theoretical dispersion curve for whistlers [Helliwell, 1965], 

(kuc/w) 2 - 1 - v 2 /w(w - « ), shown as a solid line for the independently 

y ' pe ce 

measured experimental parameters B^. Use of plane-wave, infinite medium 
theory is justified since « 1 and r plasma^ C ^ ^ [Klozenberg, et 

al., 1965]. Thus, the wave is unambiguously identified as a whistler mode 
from measurements of dispersion and polarization. 

Pulsed currents, therefore, penetrate into the plasma at the 
characteristic group velocity of a whistler wave packet, v — 3w/3k « 

O 

2c (uft /w 2 )V 2 (w « Op)- This observation is expected to hold for 
conditions where ion cyclotron phenomena are negligible, i.e. electrode 
dimensions smaller than an ion Larmor radius and pulse lengths shorter than 
an ion cyclotron period. Current systems on MHD scales are believed to be 
convected by Alfven waves [Goertz, 1980]. Most important, however, is the 
fact that time-varying currents are transported by waves, not by streaming 
particles as may be inferred from probe theory [Langmuir, 1926]. For 
example, when the electrode is biased negatively so as to collect slow ions, 
the current front is also found to propagate at the speed of electron 
whistlers [Urrutia and Stenzel, 1989]. Thus, just like in ordinary 
conductors, switched currents are associated with wave propagation effects, 
i.e. perturbations in the particle motion, rather than by the streaming 
motion of collisionless particles which is the accepted concept of dc 
currents. It is interesting to note that dc currents can be driven by 
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continuous waves through Landau damping (FIsch and Kamey, 1981; Wong et al., 
1980) but such wave -particle resonances are not necessary in transient 
currents . 

3 . 2 Moving Electrodes 

While a single pulse from a fixed electrode gives an elementary picture 
of the current flow, we now proceed to synthesize the continuous motion and 
dc current from a sequence of delayed steps and current pulses. The method 
of linear superposition, however, is only valid when the plasma parameters 
are not perturbed by the injection/collection of current. 

The concept is schematically shown in Figure 4. Ideally, we would like 

to move the electrode across at a speed v^ on the order of the current 

penetration speed along (v » 100 km/s) but this is energetically not 

feasible (^Mv x 2 > 1 GJ). We use a different yet very simple approach based 

on the linear superposition of electromagnetic fields. As shown in Figure 

4a, the continuous motion v^ is broken up into individual adjacent Increments 

(Ax “ size of electrode) and the dc current, 1^ , is divided into a sequence 

of adjacent current pulses (At “ Ax/v^ ” transit time of electrode through 

its diameter) which are applied in succesive plasma shots. For each 

electrode position, x^, and current pulse, t^, the magnetic field 

perturbations B^(x,y,z,t) are measured and stored digitally. After linearly 

superimposing all the retarded field contributions, we claim (and 

substantiate below) that the resultant field approximates that of an 

electrode moving at a constant velocity v^ and drawing a constant current 

I The expected effect of the electrode motion is sketched in Figure 4b. 

Each current pulse propagates along B^ with the speed of a whistler wave 

packet v (center frequency z inverse pulse width) while the electrode moves 

with v across B* . The result is an oblique current path or "whistler wing" 
x o 
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In the x-z plane at an angle # given by tanf - v /v^. In the moving frame, 
the wing is stationary and the current flow appears oblique to B q . The wing 
is a propagating whistler wave packet in the stationary plasma frame . 

In order to check whether linear superposition holds under the present 

conditions, the test sketched in Figure 5a has been performed. Using two 

separate electrodes at different locations, the fields B^ and Bg of current 

pulses 1^ and Ij, respectively, have been measured separately, stored, and 

later added digitally. The result is compared with the case when both 

currents are applied simultaneously and the sum field B^ ♦ Bj Is measured 

directly. As shown In Figure 5b, for the results are identical. Although 

not shown, the same is found to hold for the other two components (B , B ) at 

x y 

different field points. Thus, It is valid to assemble the "whistler wing" 
from separately measured wavelets launched by delayed pulses at different 
electrode positions. This is, of course, theoretically expected on the basis 
of Huygen's principle [Jackson, 1962]. 

We have also checked the correct temporal interference of delayed 


fields/currents to approach steady- state conditions. For a single electrode 

at a fixed position. Figure 5c shows (i) the linear superposition of fields 
_N_ 

^ \ B. from N - 20 delayed current pulses (At ■* 0.15 /ts) , and (ii) the field 
N L — i 

from a single long current pulse which is approximately equivalent to the 
superposition of short current pulses (see top trace). Except for small 
deviations at turn- on and turn-off, the temporal superposition of short, 


delayed field/current pulses correctly approximates the quasi-dc conditions. 

Having separately confirmed the validity of spatial and temporal field 
superpositions, we now show the results of simultaneous current displacements 
in space and time which approximates a moving collector (Figure 6) . For 
simplicity, we again display only one field component choosing B z since it 
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maximizes in the center of the current channel while B^, B y assume a minimum. 

Later it will be shown that B arises from electron Hall currents associated 

z 

with a radial space charge electric field and the axial magnetic field. 

Figure 6 shows contour plots of the perturbed magnetic field B^(x,z) at 

different times for a spherical electrode (diam. d - 2.5 cm « 10 r^) moving 

at an effective speed v^ - Ax/At « 2 x 10 T cm/s across and drawing a 

constant current I . ■ 1 A, The contours of constant field are inclined at 

ac 

an angle 6 * 35° “ tan 1 (v /v ) where v “ 3 x 10 T cm/s Is the observed 

X z z 

propagation speed of current pulses along B q . The obvious wing structure 

moves together with the electrode across B o as indicated by the three 

snapshots in time. The width of the wings is considerably broader than the 

size of the electrode and has a scale length on the order of a whistler 

wavelength Ay at the characteristic frequency of the wavepacket (f - At 1 - 6 

MHz, Ay “ v^/f « 5 cm). Further spreading of the wavepacket can be expected 

since whistlers can propagate oblique to B o within a ray cone of angle $ > 

11° [Helliwell , 1965]. Thus, the dc current from a moving electrode is 

entirely different from that of a stationary electrode where the current 

should theoretically be field-aligned within the flux tube of the size of the 

electrode when the plasma is collisionless. 

The difference between particle motion and current flow is even more 

pronounced when one looks at the current flow of a negative electrode which 

emits an electron beam (100 eV, 1 A, 0.8 cm diam.). It is tethered to an 

electron collector located at 15 cm across iT . Both electrodes are 

o 

electrically floating and moving together across like a tethered satellite - 
balloon system in space. Figure 7 summarizes the observations of fields and 
particles near the beam electrode for this case which slightly differs in 
density from the one shown in Figure 6. The figure displays, in a three- 
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dimensional topographic view, vs. x, z and Indicates schematically the 

beam trajectory based on light emission measurements. The insert shows the 

measured light signal of a single 50 ns electron beam pulse at two positions 

along & (Az - 50 cm) . The observed delay time (At ■■ 80 ns) yields a 

particle velocity (v b - Az/At « 6 x 10* cm/s) consistent with the beam 

energy. Since the beam electrons move much faster than the cathode, the 

angle between the beam trajectory and the magnetic field is very small [* b - 

tan’ 1 (v /v fe ) - 1.9°]. However, the whistler wing, hence current flow, is 

inclined at a much larger angle 6 c « 18° due to the slower propagation speed 

of whistlers (v « 6 x 10 7 cm/s) than beam electrons. The current is spread 
z 

out over a much wider region than the beam. Since the beam produces no 
measurable magnetic field <B X - B - B z « 0) , it is current -neutralized by 
the background electrons. The net current flows along a wave supported by 
the background electrons. The positive electrode [located outside the figure 
(at x - -5 cm)] produces a similar wing. Since the tethered electrodes are 
electrically floating, the current must close within the plasma volume. 
Although the beam strikes the chamber wall there is no net current flowing to 
the wall. Figure 7 represents a snapshot in time; the beam and wing 
translate in time along the x-direction and remain stationary in the moving 
frame as long as the plasma parameters remain constant (Ax « ± 40 cm, At - 
±2 /is). 

3.3 The Tether 

While it is well known that a bare tether wire collects charged 
particles directly from the ambient plasma [Dobrowolny, 1978; Morrison et 
al., 1978], it is usually assumed that a dc insulation decouples the wire 
from the plasma such that the only interaction is via the end electrodes. 

This concept, based on dc circuits in stationary media, needs to be revised 
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Just like the assumption that plasma currents are confined to magnetic flux 
tubes of the size of the moving electrodes. The following laboratory 
observations show that an insulated wire can excite whistler waves over 
its entire length. Consequently, the rapidly moving tether carrying a dc 
current will act like a long wire antenna. The power loss due to radiation 
of waves is, therefore, much larger than that calculated for the end 
electrodes alone. 

An insulated straight wire (0.25 mm diam. , 1.5 m length) has been 
suspended across the entire plasma column (along x, xB^) and a current 
pulsed through it (I * 1 A, At « 20 ms). With the magnetic probe, the 
perturbed magnetic field B(z,t) is measured on axis (x — y — 0) . Figure 8a 
shows typical waveforms of B^(t) at z - 10 cm from the wire in the plasma, in 
vacuum, and the difference, AB - B p - B v> which is the field created by 
induced plasma currents. The latter is a transient field of duration 
comparable to the pulse rise time (At - 1 ms) and of magnitude comparable 

2 

to the vacuum field. The shielding currents do not diffuse resistively (V a 

M q o 8/8t) but couple to a propagating wave (V a 8 /dt ) of nearly 

constant group velocity and weak amplitude decay, as indicated in Figure 8b. 

The contour plot AB - const, vs. space and time (Figure 8b) shows the 

characteristic propagation speed, v ■ 3 x 10 cm/s, dispersive properties, 

g 

v “V (w) , and wave reflections from the end anode at z “ 60 cm, vg — -0.8 x 
g 6 

10 7 cm/s . 

In order to identify the mode excited we use, as in Figure 3, a square- 

wave modulated current and an interferometer circuit to plot orthogonal field 

components B^, B^ vs. position z (Figure 9a) which yields again a right-hand 

circularly polarized wave, (B y , B x ) a (coskz, sinkz), with dispersion w(k|) 

2 2 

closely matching that of low frequency whistlers, u « k| c O e / w p (Figure 
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9b) . The whistlers propagate away from the wire in 1 z direction as is 
demonstrated by plotting fixed phase points (kz * wt ■ 0, 180 ) vs. position 
z and reference delay time r (Figure 10a) . The amplitude decay with distance 
z (Figure 9a) is mainly due to geometric spread of the wave energy from the 
exciter wire which becomes apparent from the two-dimensional phase front 
measurement of By(y,z), shown in Figure 10b. Magnetic phase fronts are 
convex with respect to the exciter with line current along x, i.e. the wave 
normals, the rays and the energy flow are diverging. The contours of 
constant wave amplitude (dashed lines) are ellipsoidal which is 
characteristic for weakly diverging whistlers launched from magnetic antennas 
[Stenzel , 1976] . 

Thus, there is convincing evidence that time-varying currents in an 
insulated wire launch whistler waves of initial wave amplitude comparable to 
the vacuum fields around the wire. As in the case of the end electrodes, one 
could now model the motion of a tether wire carrying a dc current by a 
linear superposition of wave fields from delayed pulses at shifted wire 
positions. Based on the close analogy to the moving electrodes, the result 
would again be a "whistler wing" but of thickness in x given by the entire 
length of the tether wire. Basically, the moving tether with dc current 
acts like an antenna radiating a spectrum of whistler waves with k-values 
ranging from k > 0 to k - 2*/x. . In comparison with the end electrodes, 

the tether contribution to the effective "radiation resistance" should be 
much larger, i.e. by the ratio of tether length to electrode diameter. 

Although in the optimum configuration of an electrodynamic tether the 
wire axis /, the magnetic field B^, and the velocity vector v are mutually 
orthogonal to one another, it is likely that mechanical drag, electromagnetic 
forces, and intentional maneuvers change the wire orientation. If the wire 
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tilts backwards (j ? I v) but remains orthogonal to B , an applied current 
could excite longer wavelength, low frequency Alfvenic »odes (w <. v/i, 0 < k 
< 2 w/i). However, It could lose Its radiation resistance when the wire Is 
aligned parallel to B fl . We have Investigated this non-radiating 
configuration and found that Ohmic dissipation associated with magnetic field 
diffusion will produce an effective dc resistivity for a moving wire. 

The experimental configuration for demonstrating magnetic field diffusion 
is very simple. As sketched in Figure 11a, the insulated tether wire is 
stretched along the axis of the plasma column parallel to B & with the chamber 
wall providing a coaxial return path. A pulsed current (I « 1A, At - 20 ps) 
is applied and the only significant magnetic field component, B^(r,t), is 
measured as a function of time at different radial distances from the wire. 
Figure lib shows that with increasing distance the time variation of the 
magnetic field slows down, as expected from the vector diffusion equation, 
V 2 B*- fi^a 3B/3t. Th e field penetration into the plasma is impeded due to 
induced shielding currents, J(r,t) - V x B/m q « (^r)’ 1 3(rB*)/3r, shown in 
Figure 12. The current density is distributed like a shell around the wire, 
and, in time, it spreads radially and decays in strength. Except for the 
direction of the induced current, governed by Lenz's law, the induced 
currents are the same for the switch-on and switch-off of the wire current. 


In the former case, they impede the field penetration; in the latter, they 
slow down its rapid decay. For example, Figure 13 depicts the magnetic field 
B (r) at different times t after the switch-off of the wire current. The 

4 > 

field is produced only by plasma currents. The stored magnetic field energy 


which, before switch-off, is essentially the same as in vacuum is dissipated 
by Ohmic losses, E • J - 0 ||J 2 . No waves are excited and thus no energy is fed 
back into the tether circuit (I y - 0). The conductivity a g « 2.5 (Ocm) , 
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has been calculated from the neasured field dependence B^(r,t) and the vector 
diffusion equation, 1/r 3(rdB^/flr)/dr - B^/r 2 B^ - /i^aB^/flt, and compares 
reasonably well with the Spitzer conductivity for the present plasma 
parameters (< 7 ^ ** (Oco) 1 0 kT^ m 0.4 eV) , 

Since a time varying current in a stationary wire (i Q B q ) is subject 
to eddy current losses, a dc current in a moving wire (v l B q ) will encounter 
similar losses. The field distribution may form a "diffusion wing" with 
upstream fields compressed by antiparallel plasma currents and downstream 
fields expanded by parallel currents. The Ohmic losses in the plasma either 
decrease the kinetic energy of the moving tether or are provided by the 
electrical energy source driving the tether current. 

Another important observation is the propagation speed of switched 
currents along the field- aligned wire. It is observed to be essentially the 
same as that in vacuum, v < c. Conventional transmission line theory 
predicts that the propagation speed (v a TEC) is greatly reduced by the high 
plasma dielectric constant, - 1 - Wp(w 2 - O*) " “p^c’ to v “ c/100. 
However, the shielding effect prevents the magnetic field B^ to penetrate 
rapidly for the short duration of axial current propagation [Ar - (t/^o) 1 /* 
cx 5 mm for t “ 2.5 m/c - 8 ns]. Thus, the energy flow H occurs along the 
wire surface rather than through the bulk of the plasma, which would have 
slowed down the propagation speed to that of whistlers as In the case of 
plasma currents to electrodes (Figures 2, 3). 

The two tether orientations investigated (i | i l B o > showed 
fundamental differences in plasma response (diffusion vs. wave 
propagation) . For arbitrary wire orientation or complicated structures 
(e.g. large space stations) both radiative and dissipative losses can be 
expected but may be more difficult to distinguish than in a controlled 
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laboratory environment where geometry, end electrodes, and insulation can be 
separately investigated. 

4. Nonlinear Currents 

While the study of linear fields and currents forms the basis for 
understanding electrodynamic tethers, the actual application of long tethers 
is intended to operate in a highly nonlinear regime, defined by electrode 
potentials e^/kT » 1 (typically 10 4 ) and maximum possible current densities, 
J :> ne(kT /2»rm )V 8 t Although strong currents are known to drive various 
instabilities, it is still difficult to make theoretical predictions for a 
complicated dynamic current system. The following laboratory observations 
may indicate nonlinear processes which could also occur at the contact points 
of tethers in space. 

4 . 1 Current -Voltage Characteristics 

Assuming an unlimited current supply at the negative electrode, we first 

explore the limits of electron collection by a positively biased electrode 

-¥ 

consisting of a 2 cm diam. plane tantalum disc with surface normal along B q 

and one side insulated with a mica sheet. Figure 14a depicts the time 

waveforms of the collected current at different applied generator voltages, 

V , switched on at t - 0. Figure 14b lists the voltage levels, V , 
g 6 

corresponding to the different traces a-i, and displays the voltage, 

(solid lines), measured directly at the collector. Collector and generator 

voltages differ by the potential drop along the transmission line as shown in 

the inserted circuit diagram. The voltages are referenced with respect to 

the chamber wall which forms the large negative ion collecting electrode 

(subsequently replaced by an electron beam source) . While at small applied 

voltages (V < V , » 5 V) current and voltages remain constant after switch- 
g ‘ P* 

on, both waveforms become highly time -dependent for currents on the order of 
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the Langmuir saturation current, 1^ “ Ane(kT e /2irm^) / a 4.2 A (collector 

area A - 0.78 cm*, n - 10 13 cm'*, kT - 4 eV) . For short durations (At - 1 

/is) f the current far exceeds the Langmuir limit, subsequently collapses and 

fluctuates around The large time -varying current generates a significant 

inductive voltage, Ldl c /dt, along the axial wire between generator and 

electrode (L • 3 /iH) . During the current rise, the collector voltage, V c - 

V - Ldl/dt, is reduced while the opposite occurs during the current drop, 
g 

In Figure 14c the current voltage characteristics, I c vs. V^, for the 
collector are constructed at different times, At, after switch-on. These 
curves evolve in time from an initial straight line during the current rise 
to characteristics with negative differential resistance during the current 
drop, to curves which appear like classical plane Langmuir probe 
characteristics [Chen, 1965] although they are not Invariant in time. When 
the current Is observed on a longer time scale (as shown in Figure 15) one 
finds repeated current spikes. These are not transient responses to the 
large initial current overshoot but also occur after a slow turn-on of the 
voltage which avoids the first overshoot (see Urrutia, [1987]). For repeated 
afterglows, the multiple current peaks fluctuate in timing and amplitude such 
that ensemble - averaged current traces appear much smoother than the single 
event data of Figure 15. The fluctuations may exhibit transition to chaos. 
Such study is under way although it Is outside the scope of this report. 

The observations raise two major questions, (i) how is it possible to 
collect current far in excess of the Langmuir limit and (ii) , how do the 
repeated current fluctuations arise? In order to clarify these points we 
have performed detailed in- situ measurements of plasma parameters and current 
flows to be presented below. However, prior to discussion of these results, 
we complete the current analysis by showing 1-V characteristics of an 


- 18 - 


electron-emitting electrode and a floating, tethered emitter-collector 
circuit. Figure 15 summarizes the results of different electrode 
arrangements. A 2.5 cm diam. oxide -coated plane cathode is used as an 
electron emitter. When biased to « - 200V f it emits 1^ > 10 A depending 

on plasma density (n > 5 x 10 11 cm"*). The current Is constant in time when 

e 

the cathode Is biased with respect to the large chamber wall but exhibits an 
overshoot and fluctuations when connected to a small disc electrode (Figure 
15c) . The instabilities are essentially similar to those observed when the 
disc electrode is biased with respect to the chamber wall as discussed above. 
Thus, the current fluctuations arise from the electron collection at the 
small disc electrode. In the tethered configuration, the system "floats" 
with respect to the plasma potential since no net current is drawn from the 
plasma. The electrodes assume self-consistent potentials determined by their 
Individual I-V characteristics and the applied potential difference. For 
example, during the initial current overshoot the emitter is highly negative 
(^ e « -200 V) while the collector is positive « + 50V) with respect to 
the plasma potential (^p e “ 0V) . However, as the collector current 
collapses, the collector potential becomes highly positive (^ c « + 280V) 
while the emitter potential almost vanishes so as to yield a consistently 
small emission current. As will be shown below, this behavior is the result 
of ion expulsion in front of the collector. Periodic plasma backflow causes 
repeated current peaks during which the plasma attempts to establish the 
initial unperturbed conditions. In general, the electrode which limits the 
current in the circuit charges to the highest potential. 

4.2 Plasma Currents 

In order to understand the nonlinear current system we have carefully 
measured the perturbed magnetic field, B(r,t), over the region of interest 
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and calculated the complete vector field of the total current density, 

-f 4 <4 4 

J(r,t) - V x B//i o - In general, all vector components of J are present: a 

field* aligned current, J , maximizing in each current channel, an azimuthal 

z 

4 4 

current, J., surrounding each current channel caused by electron E x B 
drifts, and a weak radial cross -field current, J r , into/out of the current 
channels. Figure 16 summarizes the spatial variation by displaying J(x,y,z) 
as vector fields in an x-y plane (Figure 16a), as field lines J projected on 
an x-y plane or throughout the volume (Figure 16b), as selected flux tubes 
within which I - constant (Figure 16c) , and as surfaces of constant magnitude 
| j*| - constant in space (Figure 16d) . Such data reveal that the current 
flows along conical spirals whose pitch length varies radially as in flux 
ropes. The radial expansion of the current channel with axial distance from 

-4 

the electrodes implies cross-field transport and an axial decrease In |J|. 

The temporal variation of the current system is indicated in Figure 17 

2 

which displays contours of constant magnetic field energy density |B| /2/i Q in 
an x-z plane at different times t after energizing the circuit. As in the 
case of linear currents (Figure 2), the field perturbation propagates along 

at the speed of a whistler wave, irrespective of the strength of the current 
or the particle speed. Although the current channel of the beam penetrates 
slightly faster than that of the background electrons, the speed of the beam 
electrons measured optically by time -of- flight exceeds that of the current 
system by more than an order of magnitude (v^ » 5 x 10* cm/s). These 
observations imply that the beam ahead of the current is current and space 
charge neutralized by background electrons as discussed earlier for small 
amplitude currents (Figure 7) . 

The current closure is a topic of general interest in current systems. 
Because the total current density, which consists of conduction and 
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displacement current, is derived via V x B//i o “ J c + dD/flt, it is always 

divergence - free . Hence the total current is always closed, whereas the 

conduction current alone may not be closed due to charging effects (7 • — 

-dp/d t, p - V • I?) . First, we note that for a substantial time (At • 200 

cm/5 x 10 T cm/s - 4 /is) the current system does not reach the end of the 

chamber, hence current closure occurs within the plasma volume free of 

boundary effects just like in space. One may realize that the fast electron 

beam rapidly reaches the chamber wall (At - 200 cm/5 x 10* cm/s - 0.4 /is) but 

since it carries no net current it has no effect on the internal current 

system. The net particle flux at the floating chamber wall must remain 

balanced (I. + I + I. - 0) , thus the plasma electron flux (I ) decreases and 
b 6 I e 

the ion flux (1^) slightly increases when the beam (1^) is absorbed at the 
wall. This is accomplished by a self-consistent plasma potential increase 
(A < kT e /e) with respect to the wall. The wall position at z - 200 cm or 
even at z • • has no effect on the current closure or disruptive 
instabilities since the current front is at z " 50 cm when they start at t 

2 /is. 

We return to the problem of current collection and closure and first show 

that the observations are inconsistent with classical dc probe theories 

assuming random particle motions and electric fields confined to sheaths. 

For example, at a collision rate i/ el - 3 x 10* s* 1 and electron Larmor radius 

r - 1.5 mm (kT - 1.5 eV, B - 30 G) , electrons can randomly cross the 
Le e o 

magnetic field in At - 1.5 /is only over distances Ar - (i/ fti At) / 2 r^ « 

3 mm, thus the maximum number of electrons which can be collected by an 
electrode of radius r would be those in a flux tube of radius r + Ar and 
length i, N - w(r + Ar) 2 in g » 5.3 x 10 1S (r - 1 cm, f - 50 cm, n g « 2 x 10 11 
cm* 3 ). The actually collected number, N - Q/e - Jidt/e » 10 A x 1.5 /is/e “ 
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9 x 10* S , exceeds those classically available, which leads to Inconceivable 
space charge problems (E^ • 2.4 x 10 # V/cm) . Thus, a more realistic picture 
is that the collected electrons must come from a large plasma volume. To do 
so, they must be accelerated across and along by Internal electric fields 
In addition, the ions should also contribute to the cross-field currents. 

The observed current closure by cross -field currents is demonstrated in 
Figure 18. Two ficticious rectangular boxes define control surfaces which 
enclose the current channels of the electron emitter and collector, 
respectively (Figure 18a) . The total current I - J j* • da* is integrated for 


each surface A and displayed vs. time in Figure 18b. One can see that the 
parallel current near the electrodes I|(z ” 5 cm) is almost entirely 
accounted for by the perpendicular current 1^ flowing through the four side 
surfaces of each box, and the field-aligned current at large distances T y (z 
50 cm) is essentially negligible for t < 2 ns. Although the radial 
current density is small compared to J| and J^, the large transverse 
surface area allows the total current to flow almost exclusively across B* . 
It is also interesting to see that the cross-field current does not 
exclusively flow through the plane between the two flux tubes but spreads 
radially in all directions. Thus, the current closure does not occur along 
the shortest path, as expected from dc concepts, but along two diverging 
overlapping whistler waves. Wave and space charge electric fields aided by 
plasma turbulence are thought to provide for the strong cross -field currents 
Unfortunately, the evolution of the strong currents cannot be followed 
for longer times, since, as noted above, another overriding effect takes 
place, i.e. an instability at the collecting electrode which disrupts the 
continuation of large current flows. This process is analyzed next. 
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4 . 3 Plasma Modifications 

The current overshoot upon switch-on of a positive voltage step to a 
Langmuir probe was observed a long time ago [Bills et al., 1962] and 
plausibly explained by the expulsion of excess ions from the sheath which 
changes from ion-rich to electron-rich [Chen, 1965]. The present 
observations extend and significantly modify this model as demonstrated by 
Figure 19 which displays the collected current vs. tine for the same step 
voltage (V - +100 Volt) to three plane electrodes of very different areas 
(3 mm 2 , 3 cm 2 , 20 cm 2 ). Since the sheath thickness (10 - 0.4 mm) is small 

compared to the probe dimensions, the duration of the overshoot (At * f pi 
0.1 ft s) should not depend on probe size if only the ions in the sheath would 
be expelled. However, the data clearly show that both the time length of the 
overshoot and the repetition time of subsequent overshoots increase with 
collector size (t a r c ) suggesting that the ion perturbation and associated 
electric fields are not confined to the sheath but to a volume determined by 


the collector size. 


Using the smallest (r • 1 mm) and least perturbing electrode at low 
voltages (V^ as < ^pi asma ) as a diagnostic Langmuir probe by the procedure 
described in Section 2, we have measured the plasma parameters in space and 


time around a large collector (2 cm diam.) exhibiting current overshoots. 
Figure 20 shows the plasma density vs. radius and time after switch-on of 
the collector voltage (V - + 100V) at Az - 1 cm in front of the collector. 
The initially uniform density profile develops a deep density depletion (6n/n 
“ 80%) inside the current channel. Plasma is expelled to the sides with 
supersonic velocities (Ar/At ■> 2 cm/2ps »■ 10 cm/s > c g ■» (kT^/m^) / 2 » 2 x 
10 6 cm/s), partly diffuses back, and is quasi-periodically re-expelled (at t 
<* 5 and 10 fts, on the average). Due to necessary ensemble averaging, the 
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repeated density perturbations are reduced in size since their timing 

fluctuates from shot to shot. Axially, the density depression is deepest in 

front of the electrode ($n/n * 100%) and gradually decreases with distance 

along (5n/n ■> 50% at z * 10 cm), implying that the unmagnetized ions (r^ 

> 20 cm) are expelled mainly across out of the positively charged current 

channel. Ions are expelled out of the region but not lost to the wall while 

electrons are attracted into the current channel and absorbed at the 

collector. In spite of the different particle flows, the space charge 

neutrality is preserved to a high degree [ (n^ - n e )/n - V • D/ne » 10 4 ], 

thus the expulsion of ions causes a proportionate loss of collected 

electrons, i.e. a current drop. After current collapse, the plasma pressure 

causes a reversal of the ion motion until the next expulsion takes place. 

The width of the overshoot is determined by the collector radius r and ion 

c 

x /2 

expulsion velocity, r /(«V c /m^) , the repetition time depends on 

r c /(kT e /m^) V 2 [Urrutia, 1987] Such dynamical effects are best observed in 
an afterglow plasma since in an active discharge they would be obscured by 
the continuous plasma production. 

Thus, in contrast to the conventional picture, the current overshoot is 
not due to an excess ion density near the collector which is removed to 
result in a uniform plasma, but to unforeseen large currents flowing in an 
Initially uniform plasma that is subsequently perturbed. This implies that 
the low time -averaged dc current cannot be explained by conventional Langmuir 
probe theories which assume an unperturbed plasma beyond the sheath. 

Likewise, the high initial current can only arise due to the penetration of 
electric fields beyond the sheath into the plasma volume. Such fields are 
measured and displayed In Figure 21. The particles experience predominantly 
the field E - caused by plasma potential gradients while pressure 
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gradients (Vp/ne) causing diamagnetic electron drifts are negligible at early 
times (t < 2 /is) and do not significantly contribute to the collection of 
electrons. Near the collector (* •• r •» 1 cm) one observes strong parallel 
fields (E * 6 V/cm) in addition to perpendicular fields (E ■ 4 V/cm) which 

Z A 

maximize near the edge of the current channel (r « 2 cm) up to large axial 

distances (z 5 30 cm) . These fields are responsible for accelerating 

electrons to velocities far in excess of the initial thermal velocity, and to 

2 

collect them over a long cylindrical surface of area (2irrz “ 380 cm ) much 

2 2 

larger than that presented by the collector's sheath <« c m 3.14 cm ). The 

ion expulsion also adds to the observed cross-field current (1^ “ 2»rrzn^e<v^> 

“ 2.4A, <v j> “ 2 x 10 6 cm/s). The large current density near the collector 

(I/wr c 2 “ nev d “ 3A/cm 2 ; n “ 2 x 10 11 cm**, v d * 10* cm/s) is associated with 

1 2 

a high effective electron temperature (^nKv > - kT >4 eV) . The parallel 

1 e eff 

electric field cannot be explained by classical collisions [ ct y - J|/E|| ” 
3A/cm 2 /(3 V/cm) “ 1 (Ocm) 1 at t “ 1 /is, a — 153 T e */ a /in A “ 123 (Ocm) 1 » 
Oj|, Spitzer, 1962]. Both the parallel drift and the x B q electron drift 
(v^ - E r /B Q “ 1-3 x 10 cm/s) far exceed the sound speed [c g - (kT^/m^) / a “ 

2 x 10* cm/s] such that ion acoustic waves are destabilized. Evidence for 
plasma turbulence is shown below. 

As the overshoot current collapses, the electrode potential, and hence 

parallel electric fields, is greatly enhanced by the inductive voltage Ldl/dt 

(see Figure 14b) . When the potential drop In the plasma exceeds the 

Ionization potential (eV^ « 15.4 eV in Ar) , and an adequate number of 

neutrals is present (n > 2.5 x 10 1S cm 8 , p > 7 x 10 4 Torr) , we observe 

o ” o 

the ionization of the neutrals by electron impact. Figure 22 shows current 
waveforms l(t) at gas pressures below (Figure 22a, p Q - 4 x 10 4 Torr) and 
above (Figure 22b, p Q “ 7 x 10* 4 Torr) onset of ionization for an applied 
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collector voltage, V c - + 100V. Ionization phenomena are usually accompanied 
by excitation of light, also displayed In Figure 22, which is easily detected 
since the surrounding afterglow plasma is dark. The ionization starts after 
the first current overshoot and produces a growing plasma cloud ("contactor") 
in front of the collector which has independently been mapped with probes 
[Urrutia, 1987]. As long as the ion production by ionization exceeds the 
rate of Ion expulsion by the electric field, the large electron current to 
the collector can be maintained. However, at t £ 25 /is both ionization and 
current decay In spite of a constant collector voltage. The finite duration 
of ionization is explained by a decrease in neutral density (n n - 3 x 10 13 
cm' 3 ) at the contactor since the fast Ion outflow (Nj-IAt/e “ 5 x 10 14 ions) 
exceeds the slow neutral inflow (N n “ 3 x 10 14 atoms for a 10 cm 3 contactor 
volume). It Is Interesting to note that with a longer repetition time (t r * 
100 /is) determined by the neutral inflow, the Ionization process repeats like 
a relaxation oscillation. 

The counters tr earning electron and ion saturation currents near the 
collector biased at » kT g /e provide the necessary conditions for the 
formation of a potential double layer [Block, 1978] . Anodic double layers 
have been observed In various laboratory plasmas [Torven and Andersson, 

1979; Alport et al., 1986] and have been investigated in the present device 
in connection with current sheet disruptions during magnetic reconnection 
[Stenzel et al., 1983]. In the presence of ionization, the double layer 
potential drop is on the order of the ionization potential. 

4.4 Microinstabilities 

Using rf probes (coax*fed cylindrical Langmuir probes, 0.1 mm diam. , 1 mm 
length) , the fluctuations in the electron saturation current have been 
measured in front (r - 0, Az •> 5 mm) of the collector and are displayed in 
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Figure 23 together with the collector current vs. time. Each current pulse 

produces a significant enhancement in the noise level demonstrating the 

existence of current -driven microinstabilities. Spectrum analysis of the 

noise [Urrutia, 1987] reveals that the instability starts in the ion acoustic 

regime (« < u> - 2* x 15 MHz) , but rapidly broadens in time to frequencies f > 
'pi 

4 f - 60 MHz, indicative of a strong two-stream instability when the 
Pi 

electron drift approaches the thermal speed. Conditions for the classical 
Buneman [1959] instability (v d > v g ) are not satisfied since the electron 
temperature apparently increases with drift velocity. 

The dispersion w vs. k of the current -driven fluctuations has been 
measured by frequency- resolved two-probe cross-correlation techniques. 

Figure 24 (insert) displays typical cross- spectral functions Cfin^z.r) x 
6 n 2 (z+Az,r+Ar)> at f - 26 MHz vs. displacement along (Az) and across (Ar) B^. 
The data indicate short correlation lengths, i.e. strong turbulence. The 


dispersion curve, obtained from correlation wavelengths at different 
frequencies, differs from the ion acoustic wave dispersion by a higher phase 
velocity (w/k| “ 3 c g ) and by the lack of a cutoff at the ion plasma 
frequency. The direction of wave propagation is toward the collector, i.e. 
with the electron drift, as expected for a current -driven instability. 


The theory of ion sound resistivity [Wesson and Sykes, 1973] predicts 
that for v “ v , kT “10 kT. and n “ 2 x 10 11 cm*’ the effective collision 
frequency for electrons is u* m 4 x 10 8 * 100 MHz, and the normalized 

wave energy density W/nkT ft “ ( 6 n/n ) 2 “ 2.5 x 10* 8 . The former value yields a 
turbulent conductivity a * » n^e 2 /m»/* « 0.56 (Ocm) which is close to the 
observed value 0 | ■» J|/E| “ 1 (Ocm) 1 at t “ 1 (is, and the latter gives a 
normalized fluctuation level 6 n/n « 5% consistent with observations. 


Thus, a possible scenario for the current overshoot is the growth of 
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current -driven microinstabilities, generating parallel dc electric fields and 
increasing the cross -field transport, which allows the collection of more 
electrons, providing positive feedback for growth of turbulence, currents and 
fields. The growth saturates when the ions begin to move out of the current 
channel in response to spacecharge and resistive electric fields. The 
strongest turbulence is observed during the current collapse when inductive 
electric fields try to maintain the current density in spite of a decreasing 
particle density, implying a high drift velocity. Electrostatic fluctuations 
near the electron plasma frequency (f - 2-4 GHz) are detected, and even 
electromagnetic radiation is detected remotely and identified as harmonics of 
the sheath-plasma resonance driven unstable by electron transit time effects 
[Stenzel, 1988] . 

The microinstabilities excited by the electron beam have been 
investigated in related earlier experiments [Whelan and Stenzel, 1985]. The 
predominant modes are Langmuir waves (w :> w^) driven unstable by the cold 
beam-plasma instability [Briggs, 1964], and re turn -current driven ion 
acoustic instabilities. The high frequency instability develops strong 
Langmuir turbulence within a distance Az “ 10 A •« 10 v^/f “ 2 cm from the 
cathode which scatters both beam and background electrons. Suprathermal 
electron tails are formed and bulk electron heating occurs [Needelman and 
Stenzel, 1987]. The scattering of Langmuir waves (w e ) off ion acoustic waves 
(u>i) produces weak electromagnetic waves :> u^, k em - k g - 

k t “ 0). None of these instabilities seriously impede the current flow. 

Since the current channel charges negatively, it attracts ions and prevents 
current disruptions by density depletions. 

5. Conclusions 

The present laboratory experiment has addressed several important issues 
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of the physics of electrodynamic tethers, l.e. tine dependence and motion of 
currents, the global current system Including propagation, closure, 
electrodes and transmission lines, nonlinear phenomena such as current 
disruptions by density depletion, anomalous electric fields leading to 
unlimited current collection, and anode contactor generation by self- 
ionlzation. It Is interesting to speculate on the implications of these 
findings on the behavior of tethers in space. 

In the linear regime, the first finding that currents from/to small 
electrodes advance at whistler wave speeds rather than Alfven speeds Is 
directly applicable to space experiments. For typical shuttle speeds (v < 

8 km/s) and effective collector sizes of D - 1.6 m diameter, the duration of 
the current pulse drawn on the collector flux tube is At - D/v — 1.25 ms - 
1/800 Hz. This time is too short to excite Alfven waves (f ^ “ 30 Hz for 0 + ) 


-s, 


but appropriate to couple to low frequency whistlers (f/f cg “10 ) . Since 

the whistler group velocity v “ 2c (uO /w 2 ) 1 / 2 “ 6000 km/s is about three 
orders of magnitude higher than the spacecraft velocity, the angle of the 
"whistler wings" will be entirely negligible (0 “0). However, this does not 
imply that the current remains field- aligned because the contact, which is 
small compared to the typical whistler wavelength (A| - v^/f “ 3 km) , 
radiates like a point source. The current spreads within the whistler ray 
cone ($ c = 19°) which results in overlapping of wings within a distance d > 

30 km along B^, possibly defining the scale length for the effective current 
closure. 

Vhile the insulated tether wire does not give rise to dc plasma currents 
it will induce ac current pulses on each flux tube. The tether should 
generate a co -moving, wedge -like whistler radiation pattern along its entire 
length (i “ 20 km). An observer on the ground may be able to detect a 


- 29 - 



whistler tone burst as the tether passes by, while a free -flying diagnostic 
subsatellite nay be able to nap the dc nagnetic field perturbations with a 
nagne tone ter. 

Electron beans injected along B^ nay be distinguished from the current 
system by the small particle spread across B^, and the higher propagation 
speed (" 40,000 km/s) which should leave the beam front current- 
neutralized. Modulating the electron beam in order to generate VLF magnetic 
fields with a large (f » 30 km) "phantom" loop does not appear feasible 
with a current -neutralized beam. It is the whistler wing structure which 
represents effectively the magnetic source distribution for exciting VLF 
waves (w “ 0) . 

A long, field- aligned, current-carrying wire moving rapidly across B q is 
subject to magnetic shielding by induced field-aligned plasma currents. The 
diffusivity (/wr)* 1 « x 2 /t and the motion v - x/t determine the scale length 
x ac (v/i^o)* 1 < 1 m for the field penetration. Hence, a long wire or a large 
wire loop may not be as effective a VLF radiator In plasmas as In free 
space . 

The nonlinear phenomena are mainly related to nicroinstabilities and the 
Ion dynamics near the positive electrode. In the absence of contactors, the 
ion expulsion by the large tether potentials (kV) far exceeds the ion ram 
energy (* 4 eV for 0 + ) , hence current disruptions and fluctuations are 
likely to occur even for a moving passive collector. However, by supplying 
the expelled ions from a plasma source, the collector can remain in 
continuous contact with the ambient plasma and draw large dc currents. The 
contactor may consist of an independent plasma source (hollow cathode) or a 
neutral gas cloud ionized by ambient electrons accelerated to the ionization 
potential in a distributed or localized electric field (double layer). In 
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dense contactor plasmas (u^/0 c » 1) the collection of electrons across B & 

Is found to be much easier than expected on the basis of classical 
collisions /® c ° n T ft )• The present laboratory 

experiment does not, however, address the effects of relative motion of a 
gaseous contactor with respect to the ambient magnetoplasma which Include 
drag, electrodynamics, plasma expansion and instabilities. 
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FIGURE CAPTIONS 


Figure 1 Schematic arrangement of the plasma device, the current system, and 
the main diagnostic tools. Two movable electrodes are biased with respect 
to each other or the wall in order to establish a pulsed current flow 
through a source -free afterglow plasma. 

Figure 2 Perturbed magnetic field, B^, vs. time, t, at different axial 
distances, z, from an electrode drawing a pulsed electron current I^(t) 
(bottom trace) . The magnetic field develops transient oscillations whose 
frequency decreases with arrival time, characteristic of low frequency 
whistler wave packets. 


Figure 3 Square -wave modulation (f - 1 MHz) of the collector current 
generates cw whistler waves which are analyzed In an Interferometer circuit. 

(a) Perturbed field B x vs. axial position z (|B^) for different delay times r 
In the signal path. The data show undamped waves at f - 1 MHz propagating in 
+z direction with phase velocity v ^ « 2 x 10 ? cm/s (wavelength A| « 20 cm). 

(b) Orthogonal magnetic field components B (z), B (z) exhibit a 90° phase 

x y 

shift in space which indicates right-hand circular polarization (demonstrated 

in insert) . (c) Normalized dispersion relation, experimentally obtained from 

wavelength measurements at different frequencies (dots) , and theoretically 

given by k|| 2 c 2 /w 2 - 1 - w*/w(w- O g ) for the given experimental parameters 

f - 28 MHz, f « 8050 MHz (solid line). Both polarization and dispersion 
ce pe 

measurements demonstrate that current perturbations excite whistlers. 


Figure 4 Schematic diagram demonstrating the measurement method for 



obtaining currents from moving electrodes, (a) The continuous motion of an 
electrode is replaced by a sequence of displacements, Ax, while the 
continuous current, I. , is broken up into a sequence of pulses delayed by At 

GC 

consistent with the velocity v - Ax/At. (b) As the electrode moves with 
velocity v^ across it emits at each position a current pulse which 
propagates along at the whistler wave speed. The superposition of all 
pulses yields a wing- like current pattern at a Cerenkov- type 0 - tan 1 
(v /v^) . In the experiment the field distributions, B^(x,z), are recorded 
from repeated shots for i * 1...20 different electrode positions and 
subsequently superimposed. 

Figure 5 Experimental tests to establish the validity of linear field 

superposition required for constructing the whistler wings, (a) Schematic 

arrangement of two electrodes pulsed either separately or simultaneously and 

the resultant magnetic field perturbations, B^. (b) Measured waveforms of 

the current pulses, I(t), and the magnetic field components, B z (t), for each 

electrode separately (middle traces), which, when added (bottom trace), 

and found to be identical to the field generated by applying simultaneous 

current pulses to both electrodes, (c) 20 delayed, short (0.15 ps) 

current pulses added in order to simulate a long pulse (which could be 

extended to approach I. ). The measured fields, B (t) of the short pulses 

ac z 

are also added and superimpose closely to the field of the long pulse (quasi - 
dc) indicating proper measurement of the Fourier components of the wave 
packets . 

Figure 6 Contour plots of the perturbed magnetic field, B (x,z) (arbitrary 

z 

units, equal contour spacing, measured for x > 0, z > 5 cm), for three 
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different time steps or positions of the electrode moving at a constant 
effective speed v^ across B^. The electrode drawing a constant current (I^ c “ 
1 A) generates a co*moving "whistler wing." The currents are in steady-state 
in the moving frame of the electrode, but in the rest frame of the 
plasma the currents are transient wave currents. Whistler wings are also 
generated when the electrode collects ions or emits electrons. When scaled 
to tethered electrodes in space, the whistler wing is expected to be field- 
aligned (v « v w ) , but the wave/current spreads from the point 
radiator/electrode within the ray cone (0^ « 19°) rather than forming a long 
field-aligned path into the ionosphere. 

Figure 7 Current and particle wing of an electron beam source moving with 

constant velocity v^ «* 200 km/s across B o - 20 G. Contours of constant 

magnetic field perturbations, B z (x,z) (arb. units, equal contour spacing), are 
shown at a fixed time when the beam source is at x “ 10 cm (the co-moving 
positive electrode is displaced by Ax - 15 cm, forming a wing as in Figure 6). 
The whistler wing which carries the net current (dashed line) is clearly 
separated from the particle trajectory (dotted line) which is essentially 
field-aligned (v^ » v ) . The speed and presence of the beam particles are 
measured optically by time -of -flight. The insert shows the light emission 
from a 50 ns, 100 eV, 1A beam pulse at two positions separated by Az - 50 

cm. The measurements demonstrate that the beam is current-neutralized by 

plasma electrons which transport the current as a collective perturbation, 
the whistler mode. Consequently, a modulated beam cannot be used as a long 
magnetic antenna. 

Figure 8 Magnetic perturbations created by a 1 A, 20 ps current pulse 
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through an insulated wire stretched along x (±B o > through the plasma, (a) 

Perturbed magnetic field B y (t) at a distance z - 10 cm from the wire 

measured in the plasma (B ) , in vacuum (B ) , and the difference (AB - B - 

yp yv y jr 

B ), created by induced plasma currents, (b) The field perturbation, ABy(t) , 
at different axial positions exhibits propagation delay, weak damping, and a 
reflection at z «• 60 cm where the end anode is located. Thus, the 
perturbation propagates as a normal mode of the plasma. 

Figure 9 (a) Interferometer traces of two orthogonal field components, B^(z) , 

By(z), vs. axial position, z(|B o ), for an ac current (f - 4 MHz) in the wire. 
The 90° phase delay indicates right hand-circular polarization (see Figure 
3b). (b) Normalized dispersion relation for parallel whistler waves (solid 
curve) and measured dispersion (dots) , showing that the magnetic 
perturbation is a whistler mode. 


Figure 10 (a) Wave maxima (filled circles) and minima (open circles) versus 
axial position, z, at different delay times, r, in the interferometer 
reference path, showing that whistlers are emitted with phase velocity V ph “ 
dz/dr - f A|j « 4 x 10 7 cm/s in both directions + z along B^ from the wire at 
z - 0. (b) Two-dimensional phase fronts (open circles - minima, filled 

circles - maxima) and contours of constant relative amplitude (dotted lines) 
of whistlers emitted by the wire at y - z - 0. If the wire would move with 
velocity v^ i and carry a dc current it would generate a whistler wing 
similar to the end electrodes but along its entire length, Consequently, 

a moving dc circuit is subject to radiation losses. A fixed observer should 
detect a whistler tone as an electrodynamic tether passes by. 
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Figure 11 (a) Experimental arrangement to investigate the interaction of an 
insulated, field-aligned wire carrying a pulsed current, I(t), with the 
ambient plasma, (b) Time variation of the azimuthal magnetic field, B^(t), at 
different radial distances r from the wire. The magnetic field diffuses 
into the plasma rather than exciting a normal mode of the system such as a 
magnetosonic wave. 

Figure 12 Radial profile of the field-aligned induced plasma currents, J z - 

V x B /n , which impede the rapid penetration of magnetic fields into the 

plasma. With time t after switch-on/off the coaxial current shell (see 

inserted schematic) increases in radius, broadens in width and decays in 

amplitude. The induced electric field has the same distribution as the 

current density and is obtained from Ohm's law, “ oj E^, with the 

2 -* 

conductivity <T| determined from the vector diffusion equation V B ” 

Figure 13 Magnetic field, B^(r), at different times, t, after switch-off of 
the wire current. The initially stored magnetic field energy is entirely 
dissipated in the plasma by Joule heating. A dc current-carrying wire moving 
at a velocity v y l B q would induce shielding currents (antiparallel to I dc on 
upstream side, parallel on downstream side) which continuously dissipate 
energy of the primary current circuit. 

Figure 14 Nonlinear current-voltage characteristics of an electron collector 
including inductive circuit effects, (a) Collector current, I^(t) , vs. 
time, t, at different applied pulse generator voltages, V g [see scale in 
(b) ] . (b) Collector voltage, V c (t), vs. time, t, at different generator 
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voltages, V (t) . The difference, V - V , is mainly due to the inductive 
g c g 

voltage drop, Ldl c /dt, on the transmission line of the circuit (see schematic 

insert), (c) Current -voltage characteristics, I vs. V , at different times, 

c c 

t, constructed from the waveforms I (t) and V (t). The classical plane 

c c 

probe, thin sheath I-V characteristic should saturate at 1^ - 

i /» 

Ane(kT /2mn ) ' « 4.2 A while the observed curves exhibit time- dependent 

© © 

I-V curves first without saturation, then with negative resistance, and only 

at the temporal current minimum resembling the theoretical Langmuir probe 

characteristic. The peak current, I , scales almost linearly with 

c , max 

applied voltage, V , but eventually saturates for I c aav £ 100 A due to 

g 1 

limited ion- saturation current at the return electrode (chamber wall). 

Figure 15 Characteristics of a current system established with two floating 

electrodes emitting and collecting electrons on different flux tubes. 

(a) Schematic set up of electrodes (radius r “ 1 cm, spacing d « 8 cm l 

- 30 G) and major diagnostics to measure currents In the plasma. (b) 

Emitted electron current, I (t), exhibits disruptions and fluctuations when 

e 

connected to the electron collecting disc but remains constant when 
connected to the chamber end wall. The unstable currents are created by 
plasma modifications near the collector, not the emitter, (c) Charging of 
the "floating" electrodes in an unstable current system. During the initial 
current "overshoot" the collector is slightly above plasma potential while 
the emitter is charged highly negative (” -200V) injecting a fast electron 
beam. Upon current disruption the collector charges to high positive 
potentials (+300V) while the emitter is close to the plasma potential. 
Subsequent current peaks repeat the sequence. In general, the electrode 
which limits the current flow develops the highest charge. 
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«4 

Figure 16 Spatial characteristics of the current density, J - 7 x H, in the 
plasma with immersed electrodes at x - + 4 cm, y - z - 0. (a) Vector field 

of perpendicular (J^, J^) and parallel (J^) current density in an x-y plane 
at a fixed time t after energizing the current system (100 V, 10 A). The 
transverse currents are both due to electron E x ? drifts , caused by space 

*4 2 

charge fields, and electron diamagnetic drifts Vp x B 0 /(neB 0 ) associated 

with heating. The field-aligned currents (normal to x-y plane, but plotted 

in + y direction for purpose of display) are due to electron injection and 

collection. (b) Field lines of the current density vector j*(x,y,z) in an x-y 

plane (left) and in x-y-z space (right). The lines are nested helices of 

varying pitch lengths (pitch i -* «> on axis, i -»0 at boundary of flux 

z z 

tubes), (c) Two characteristic current "tubes" whose axes follow the J line. 

4 -4 

The cross sections conserve the current I - J -A. With distance from the 
electrodes the spiral radius and tube cross section expand due to a 
decreasing current density. Note, that particle orbits differ from current 
lines. (d) Surfaces of constant current density, |J*| - constant, at a fixed 
time in space. Such surfaces bound regions of higher current density but do 
not imply lack of current closure since currents can cross the surfaces. 

Figure 17 Temporal properties of the current system, showing contours of 

2 - 0 8 

constant energy density (AB /2 m q “ 1.2 x 10 J/cm per contour) of the 
perturbed magnetic field in the x-z plane at different times t after switch- 
on of the electrode potential V. Even for strong currents (I » 10 A, AB/B q “ 
3%) the magnetic perturbation, i.e. net current, is transported by a 
collective plasma mode (a whistler wave) and not at the speed of the beam 
electrons (v^ “ 5 x 10* cm/s) which is measured optically. The beam front 
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which produces no measurable magnetic field must be both current-neutralized 

(J - 0 since B - B "0) and space -charge neutralized (p-0, E - 0, J , - 
z * y ' 

0 since B - 0) by background electrons. 


Figure 18 Current conservation and closure analysis, (a) A ficticious box is 
drawn around each current channel and the net current, I - Jj* • da, through 
each control surface is calculated. The field aligned currents are 1^, Ig 
for the electron emitter, and I^q, 1^ for the electron collector box. The 
cross-field currents are 1^, 1^, 1^, Ig for the beam, and Ig, 1^, Ig, 1^ for 
the collector. (b) Currents I. through different surfaces for each box 

1 A- 11 

versus time prior to disruption. Current conservation “ 0. 1^ “ 0 

i-1 i-6 

is satisfied. The parallel currents near the electrodes are for t < 

2 ps mainly accounted for by perpendicular currents 1^ •> ^ 1^, I^q « Ig + 

-2. i-3 

y I^, and not by distant field-aligned currents 1^ « 1^, 1^ « I^q- The 

perpendicular currents do not only flow through the joint interface between 

_ 2 _ 

the two current channels Ig < ^ J i “ J i’ ex P ected from dc concepts to be 

i-3 i-7 -» 

the path of least resistance, but flow in all directions across B o> 

consistent with a diverging wave current. 


Figure 19 Collector currents, I(t), for three electrodes of different sizes 
biased by the same step voltage, V(t) - +100 V (Atp U i se > t r £ se « 50 ns). 

The upper three traces show the long-time behavior, the lower traces show the 
first overshoot on an expanded time scale. The results indicate that both 
the width and repetition time of current overshoots increase with collector 
size, hence cannot be explained by ion expulsions from a sheath but from a 
volume determined by the size of the collector. The current rise time is 
influenced by global circuit properties, e.g. the inductance of the feeding 
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transmission line (L * 3 /iH) . 


Figure 20 Radial density profile, n fi (r), versus time, t, after switch-on of a 
large collector voltage (V^ — + 100V) . The density Is obtained from 
complete Langmuir probe traces of a small (r * 1 mm) probe scanned at Az - 1 
cm in front of the 2 cm diam. plane collector. The initially uniform 
density profile develops a deep depletion coincident with the disruption of 
the collector current. Plasma expelled to the sides of the current channel 
flows partly back which allows repeated smaller current overshoots. 

Figure 21 Vector map of the electric field in an x-z plane at a fixed time 
during large current flows and beginning density depletions near the 
electrode. The electric field, E - -V^, is obtained from determining the 
plasma potential, at each point from complete Langmuir probe traces. The 
classical Debye sheath thickness is s » 10 * 0,4 mm, the classical Ohmic 

parallel field is E| - J||/ 0 | ** 0.024 V/cm, indicating that the observed 
large fields (E| < 6 V/cm) outside the sheath are anomalous. The 
perpendicular space charge electric field, E^ < 4 V/cm, gives rise to electron 
Hall currents (see Figure 17) and expulsion of the unmagnetized ions from the 
original current channel. 

Figure 22 Collector current and light emission from the plasma in front of 
the collector vs. time at low neutral pressures (top traces) and high 
pressures (bottom traces) at the same collector voltage (V^ - +100 V) . The 
energization of electrons by electric fields in the current channel gives 
rise to ionization at high neutral presssures, accompanied by light 
emission. Replenishing the expelled ions by ionization allows sustained 
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large current collection by this self- generated plasma "contactor." When 
the ion outflow exceeds the neutral inflow (t £ 30 /is ) the ionization and 
current collection disrupt and go over into a slow relaxation node (t f * 
100 /is) . 

Figure 23 Collector current and density fluctions vs. time showing evidence 
for current -driven microinstabilities near the collector. 

Figure 24 Normalized dispersion relation of density fluctuations in the 

current channel at Az « 1 cm in front of a 2 cm diam. collector drawing I c « 

10 A at V “ +100 V. The noise from two movable probes is analyzed with 
c 

tuned narrowband amplifiers and cross-correlated which yields wavelengths 
along and a nonoscillatory correlation curve across (see examples 
inserted) . The fluctuations propagate along B^ in the direction of the 
electron drift with phase velocity higher than the sound speed (wA| > 

3 c g ) and at frequencies higher than the ion plasma frequency. For 
comparison the theoretical ion acoustic dispersion is shown (dashed curve) , 
indicating that at large electron drifts (v d £ v^) sound turbulence 
degenerates into two-stream turbulence. 
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